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Hydrodynamic Effects of Gales on Lake Erie 

Garbis H. Keulegan 

The coefficients of wind stress and the sea roughness are derived from records of water 
levels and wind intensities relating to gales passing over Lake Erie during the past fifty years. 
The reduction of wind velocities measured at various elevations in four cities on the lake 
shore to a single "effective wind velocity" over the lake is explained. A theoretical determi- 
nation of the wind tides agrees well with the observations. The coefficient of wind stress 
derived from the observations is nearly the same as that obtained by Neumann from obser- 
vations in the Gulf of Bothnia. Both the coefficient of wind stress and the sea roughness 
are found to decrease with increasing wind velocities. 



1. Introduction 

For a wind blowing steadily with a uniform in- 
tensity over a lake, the relationship between the 
wind intensity and the frictional stress induced at 
the water surfaee may be found by considering the 
wind tides produced at the ends of the hike. Wind 
tides are those displacements of the water level from 
the initial undisturbed level that are due solely to 
the action of the wind. The relationship involves a 
characteristic quantity called the coefficient of wind 
stress. There will be an (dement of uncertainty in 
the coefficient derived if the wind stress correspond- 
ing to a wind of constant intensity is a function of 
distance, or fetch, and only the wind tides at the 
lake ends are known. In such cases the analysis 
yields average values of the wind stresses, and if it 
be true that the stress is also a function of fetch, the 
derived coefficient of wind stress will be different for 
lakes of different lengths. 

Another source of uncertainty in the values of the 
coefficient of wind stress is in the derivation of the 
formula which relates the wind stress to the wind 
tides. In working out the formula, due attention 
must be given to the nonuniformity of depth of water 
in the lake and of width of the lake surface, and also 
to the magnitude of the friction at the bottom of the 
lake. The turbulence pattern in the waters of a lake 
under the action of wind is at the present time far 
from being depicted. The partition of velocities in 
the drift current at the upper layers and in the re- 
turning gravity current in the lower layers, and the 
magnitudes of the mixing lengths, are questions which 
have not been studied. Lack of information of this 
nature prevents one from appraising accurately the 
the frictional force at the lake bottom. Besides, if 
the boundary of a lake exhibits very marked irregu- 
larities and there are also numerous islands distrib- 
uted in the regions close to the lake ends, the 
derivation of the exact wind tide formula is very 
difficult and in some cases may be impossible. These 
difficulties are present to some extent in the case of 
Lake Erie. 

A second problem in connection with the action of 
wind on lakes is the relationship between the effective 
roughness of the sea and the velocity of the wind. 
Meteorologists regard the Prandtl law of velocities 
as representing satisfactorily the state of movement 
of the wind in the layers of air above the water 



surface. This law requires a roughness length which 
is called the roughness of the sea. Once the coeffi- 
cient of wind stress based on the wind velocities at 
points of known distance from the surface of the 
water is known for a given lake, the sea roughness is 
readily obtained from the value of the coefficient. 
Since sea roughness is an inferred quantity and is a 
function of stress, the likelihood exists that the sea 
roughness is a fund ion of fetch. If the determination 
of roughness is made on the basis of average stresses 
at the lake surface, the value of roughness may vary 
for a given wind velocity for lakes of different lengths. 
It is proposed in the present work to obtain the 
values of the coefficient of wind stress and of the sea 
roughness for Lake Erie. The determinations are 
brought about by considering the maximum wind 
tides observed during the severe westerly gales pass- 
ing over Lake Erie in the past half century. 

2. Sources of Original Data 

The hydrographic charts of the United States Lake 
Survey showing the surfaee changes of Lake Erie 
during gales constitute the chief source of material 
for the present investigation. These charts show 
the water stages at the extremities of the lake, and 
the values of the hourly movement of wind and its 
direction for a few stations on the southern coast. 
Two examples of such charts are given by Freeman 

The stations for which charts of water stages are 
given have been changed from time to time. Buffalo 
and Amherstburg, or Gibraltar, are the localities 
near the ends of the lake for which the water stages 
are given for westerlies during the earlier years. 
For easterlies the corresponding localities are Toledo 
and Buffalo. In these charts complete wind data 
are shown for Toledo, Cleveland, Erie, and Buffalo. 
In the charts of later years, water stages are shown 
for Toledo and Buffalo. However, the representa- 
tion of wind data has not always been complete, 
and in such instances the missing information has 
been received from the U. S. Weather Bureau. 

It was desired to handle the analysis on the basis 
of the wind tides observed at Toledo and Buffalo. 
The data for Toledo, where missing, were deduced 
from the data for Gibraltar and Amherstburg. This 
is feasible as in general a correlation exists for the 
wind tides of nearby localities. The desired correla- 



tion factor could be determined, for in some of the 
charts the simultaneous values of the water stages at 
Toledo and Amherstburg are given. In this de- 
termination use was made also of Blunt's work [2]. 
In the determinations of the coefficients of wind 
stress, the observed wind tides must be corrected for 
barometric pressure changes. Data needed are the 
hourly values of the barometric pressures at Toledo 
and Buffalo for the periods of the gales considered. 
For these the U. S. Weather Bureau has kindly 
allowed us access to the proper records. 

3. Determination of Effective Winds 

3.1. City Wind Velocities at a Common Elevation 

The evaluation of the coefficient of wind stress will 
be made with respect to wind velocities at a height of 
25 ft above the water surface of the lake. As these 
velocities were not observed and the only available 
velocities are those measured over the cities of Toledo, 
Cleveland, Erie, and Buffalo, the lake wind velocities 
will be deduced from the city velocities. The re- 
duction to the desired values will be made in three 
steps. 

As the anemometer heights during a given period 
have been different in the cities mentioned, or have 
been changed in a given city from time to time, the 
first step of the analysis is to reduce all the measured 
wind velocities to the values appropriate to a common 
elevation. The elevation selected was 165 ft above 
the ground. 

The region of air above the ground may be con- 
ceived to consist of two layers, the boundary layer 
below and the layer of frictional influence above. 
The boundary layer is relatively thin, and the move- 
ment of air in it is controlled entirely by frictional 
stress as horizontal pressure gradient and deflecting 
forces may be neglected. In a vertical traverse the 
frictional stress is constant, and the direction of the 
stress vector coincides with the direction of the wind. 
For this layer, Prandtl [3] gives the law of velocities 



—=5.75 log-? 



(1) 



where z is the elevation of a point at which the velo- 
city is Uj u^ is the shear velocity, and e is the effective 
roughness. Denoting the frictional stress at the 
ground by r and the density of air by p a , 



U* = -\'T(ilpa. 



(2) 



When the actual roughness is large, the reference 
level from which to measure the height z is a matter 
of importance. It may be that in cities with an 
Average height of buildings of 2K the reference level 
can be taken at the distance K from the ground. 
Then, following Prandtl, the effective roughness e 2 
of a city may be taken as 



€2=2^/30. 



(3) 
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Figure 1. Relation of wind velocity, u, to height over the city of 
Tokyo, z, from measurements by Siotani and Yamamoto. 



We may now consider, to test the utility of the 
rules given above, a known case of measured veloci- 
ties over a large city. In the autumn of 1946 
Siotani and Yamamoto [4] obtained the air velocities 
over the central part of Tokyo at elevations of 60, 
55, 45, 35, and 25 m above the ground, using hot- 
wire anemometers on the tower at the Central 
Meteorological Observatory. The neighboring- 
buildings being 10 to 20 m in height, K in this case 
was 7.5 m. 

In figure 1 are reproduced some of the velocity 
values from the original curves of these authors. 
In the plot z is the height in meters measured from 
the reference level of 7.5 m above the ground. The 
straight lines drawn through the points meet at a 
common point on the abscissa axis, log e 2 =— 0.3. 
Thus the effective city roughness is e 2 =0.5 m, a 
result satisfying the rule implied by eq 3. Denoting 
by u , the value of the velocities to be read from 
the straight lines of the figure for 2=1, the shear 
velocity may be computed as %* = — (u )/(5.75 log 
e 2 ). The shear velocities thus computed are also 
entered in figure 1. Representing the data in 
figure 1 in dimensionless form, the results are the 
three points in figure 2. The distribution of the 
points confirms the relation given by eq 1, and this 
fact points to the possibility of applying the rule to 
cities in general. 

Estimating the average height of buildings in 
American cities as 30 ft, the effective roughness may 
be taken to be 1 ft. Measuring the elevation of 
points z from a reference level 15 ft above the 
ground, denoting the position of an anemometer 
by z a , the measured velocity by u a , the standard 
height by z 2 (=150 ft), the velocity at this height by 
u 2 %, the roughness of the city by e 2 (=1 ft), the 
formula to reduce the anemometer velocities to the 
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Figure 2. Relation of wind velocities measured over the city 
of Tokyo to the Prandtl law. 

velocities at standard heights is 

?=log?/logJ- (4) 

Ua €2/ e 2 

The curve in figure 3 is established through this 
relation and is used to reduce the anemometer 
velocities to the velocities corresponding to the 
standard height z 2 ( = 150 ft). For the purpose of 
comparison a second curve is shown, evaluated with 
e 2 =2 ft. The disparity between the curves being 
small, it is seen that errors arising from the uncer- 
tainties in the exact value of the city roughness 
are not very appreciable. 
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FlGUBB 3. Ratio of wind velocity at a standard height of 150 
feet above a city, u^, to the wind velocity indicated by an 
anemometer, u a , as a function of the height of the anemometer, 
z a , and the roughness of the city, e 2 . 



3.2. Wind Velocities Over the Lake Surface 

The second step is to infer the wind velocities over 
the lake surface from the wind velocities u 22 corre- 
sponding to the standard height z 2 (=150 ft) over 
the city. To establish the desired values resort is 
to be made to the characteristic relationships in- 
volving the various quantities of the layer of fric- 
tional influence. Kegarding the hydrodynamics of 
the layer, the analysis of Rossby [5] proves to be 
particularly fruitful. 

The quantities that will be of value here are the 
height, i7, of the boundary layer, the wind velocity, 
U H , at the upper surface of the boundary layer, the 
gradient wind, U g , and the angle <t> s between the 
direction of the gradient wind and the direction of 
the wind in the boundary layer. The Rossby 
relations are 



U H = U g I cos <f> s = sin </> s h 

y2 



(5) 



nud 



11=0. 54 F -^ sin fa, k 0.065, (6) 

log £/,//«= 1.694 cot (/>, log sin fa+ 1.441 . (7) 



Here /represents the quantity/=2Q sin X, 12 being the 
angular speed of the earth's rotation, and X, the lati- 
tude. The roughness of the surface below is e. The 
importance of the effective roughness, of the ground 
or of the sea, for the quantities pertaining to the lower 
layers of the atmosphere can easily be seen. 

Once it is assumed that the magnitude and the 
direction of the gradient wind are the same over the 
lake water and over the coastal regions next to the 
lake, eq 5, 6, and 7 provide the means by which a 
relationship may be established connecting the ve- 
locity of the wind over the surface of the lake with t lie 
velocity of the wind over a city on the borders of the 
lake. 

For a second definition of terms the diagram in 
figure 4 may be consulted. Assuming a gradient 
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Figure 4. Notation diagram for winds over a lake and a city 
on its shore. 
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wind constant both in intensity and in direction 
blowing in the higher regions over the lake and over 
the city, let Hi denote the height of the boundary 
layer over the lake and H 2 that over the city. Let 
€1 denote the effective sea roughness and € 2 the rough- 
ness of the city. Let u x denote the wind velocity at 
the upper surface of the boundary layer over the 
lake and u 2 that over the city. Let u 22 denote the 
wind velocity over the city at the standard height 
z 2 (=150 ft) and u n that over the lake at the standard 
height Zi (=25 ft). Let 4> s2 be the angle between the 
direction of the wind over the city, and the direction 
of the gradient wind, and <j> s i the corresponding angle 
over the lake. Then, according to the Prandtl law 
of velocities given by eq 1, the formula which reduces 
the wind velocities over the city to those over the 
lake is 



Un_ 
U 22 



~u 2 



log - / log 



log 






€l , 



^2 
€2 



*2 



(8) 



the computations being made through eq 5, 6, and 
7. The results are shown in the form of curves in 
figure 5, with the ratio Vn/n 22 plotted as function of 
u 22 for various sea roughnesses. It is seen that the 
ratio u n /u 22 for a given wind velocity depends to an 
appreciable degree on the sea roughness. In using 
the graphs for the desired reductions of wind veloci- 
ties for a given storm, it was necessary to assume a 
value for the sea roughness. The value selected was 
€! = 0.3 cm. Evaluations from the data of wind tides 
given later in this paper confirm this chosen value as 
an average. 

It will be remembered that the curves of figure 5 
are based on a city roughness of value e 2 =l ft 
(30.5 cm). Accordingly, the figures show that if the 
sea^roughness be €i = 0.3 cm, the wind velocity ob- 
served at the lake surface at a height 25 ft above the 
water will be about 1.12 times the city wind velocity 
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Figure 5. Relation between winds at a standard height of 25 
feet over a lake, Un, to winds at a standard height of 150 feet 
over a city, u 2 - 2 , for various sea roughnesses. 



observed at a point 165 (=150+15) ft above the 
ground for the higher velocities of the city wind. 
It will be found, considering the data of figure 3 also, 
that if the city wind velocity be observed at a height 
315 ft above the ground, it gives the wind velocity 
over the lake. This is in agreement with the findings 
of Kadel [6] in an investigation of wind observations 
at Buffalo. 

The roughness of the sea being much smaller than 
the roughness of a city, there is a change in the wind 
direction in passing from a city over a lake surface. 
The change is (0 s2 — </> s i), and the estimated values 
as deduced from the Kossby relations are shown in 
figure 6. The magnitude of change depends on the 
wind velocity and the sea roughness. If the sea 
roughness be chosen as ei = 0.3 cm, the value of the 
average change is about ($ s2 — <t> s \) — 7° for the higher 
wind velocities. 

3.3. Effective Lake Wind Velocity 

The third step is to determine the effective lake 
wind velocity, which is the velocity of the wind 
necessary to produce the observed wind tide under 
the condition that the wind blows with uniform in- 
tensity along the axis of the lake. When the actual 
winds at Toledo, Cleveland, Erie, and Buffalo during 
a storm are examined, it will be seen that these winds 
differ from each other in intensity and also in direc- 
tion. The need for a formula to give the effective 
wind is readily seen. 

The outlines of Lake Erie are shown in figure 7. 
The direction of a wind blowing from the south is 
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Figure 6. Difference in direction of winds over city and lake 
as function of wind velocity over the city for various values of 
the sea roughness. 
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CLEVELAND 



Figure 7. Sketch map of Lake Erie showing locations of cities 
considered and notation for directions. 



chosen as the reference direction from which to 
•measure wind directions. The lake axis is drawn so 
that it divides the area of the lake into two equal 
parts. Its direction is approximately a =68°. In 
| the cross section through the point x on the lake axis 
the velocity of the wind at an elevation of 25 ft 
above the water surface is u U} and its direction is 0, 
an angle less than 180°. Denoting the effective lake 
I wind by V and the length of the lake axis by Z, 



1 C L 



V*=- 



cos (6—a)dx. 



(9) 



This formula is applicable when (B—a) is small and 
the resulting wind velocity is large, for then the effect 
of the earth's rotation may be neglected and it may 
be imagined that the resultant movement of the 
water is in the direction of the lake axis. Putting 
,M=u n /u 22 , where u n /u 2 2 is determined by eq 8, 
eq 9 becomes 



-if 



M 2 ui 2 cos (e—a)dx. 



(10) 



In this transformation no provision is made for 
the circumstance that the wind directions are not 
the same over the lake and over a city on the shore 
of the lake. Theoretically, in the above, denoting 
the difference (0— a) as 0, the factor cos /3 must be 
replaced by cos (j3+y) where 7 is the difference 
(0s2— <t>si)- Now cos (£ + 7) = cos jS cos 7— sin ft 
sin 7, where 7 is a small quantity not exceeding 7°. 
Also the larger and significant values of wind tide 
are associated with winds for which the absolute 
values of are less than 7r/4. Therefore it is readily 
seen that the resulting error is small when 7 is 
neglected. 

When the exact variation of u u ~Mu 22 as a function 
of x is not known, the formulation must be based on 
discrete values of wind velocities. Let V T , V c , V Ej 
and V B be the wind velocities u 22 at Toledo, Cleveland, 
Erie, and Buffalo. The corresponding values of M 
and are denoted by the same subscripts. Let the 
lake axis be divided into three segments of equal 
length. It may be supposed that the winds pre- 
vailing over the end segments are deducible from V T 
and V B . The middle segment may be divided into 
two equal parts in which the winds are deducible 
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Figure 8. The reduced data for wind tides, AH, and effective 
winds, V, for storm of December SI, 1911. 

from V c and V E . Thus instead of eq 10, we may 
write 



V 2 =\[M%V 2 



\[M%V 2 T cos(d T -a) + M 2 B V 2 cos(d B -a) + 
hM%Vl cos(0 c -a) + fM|F| cos(^-a)]. 



Now for the higher wind velocities, in which we 
are interested, the variation of M is not large as 
may be seen from the curves in figure 5. Therefore 
without significant loss of accuracy, we may write 



V 2 =M 2 Vl 

V%=\[V* T cos (6 T - a) + V% 

IT/2 ™ (a «AJ_IT/s 



„ cos (6 B —a) + 
\V% cos (d c -a) + iVI cos (6 E -a)]. (11) 



In this represent *tion V 2 is, in a sense, the effective 
city wind velocity for the southern coast of the lake 
over the entire length between the lake extremities. 
Values of M as a function of V 2 are shown in table 
1 . The table is established on the basis of the curve 
in figure 5 associated with €i = 0.3 cm, after identi- 
fying u 22 as V 2 and u n as V. 



Table 1. 



Dependence of the multiplying factor M on the 
effective city wind velocity V% 



V 2 

V 2 


M2 


vl 


M2 


(mph) 2 




(mph)* 




100 


1.150 


1,300 


1.254 


200 


1.175 


1,400 


1.257 


300 


1.193 


1,500 


1.260 


400 


1.205 


1,600 


1.262 


500 


1.214 


1,700 


1.264 


600 


1.223 


1,800 


1.266 


700 


1.229 


1,900 


1.268 


800 


1.234 


2,000 


1.270 


900 


1.239 


2,100 


1.272 


1,000 


1.243 


2,200 


1.273 


1,100 


1.247 


2,300 


1.274 


1,200 


1.251 


2,400 


1.275 
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The method of velocity reduction explained above 
was applied to all the wind data of the storms under 
consideration. The final values of effective wind for 
the storm of December 31, 1911, are shown in figure 8 
as a matter of illustration. 

4. Determination of Wind Tides 

If the nodal point of the displacements of the sur- 
face water in Lake Erie during a storm be known and 
also the level of the undisturbed water prior to the 
storm, the coefficient of wind stress may be obtained 
by considering the wind tide at either end of the lake. 
Practically, however, this method of evaluation has 
to be avoided as it is difficult to determine the nodal 
point and the level of the lake with exactness. It is 
more satisfactory to base the determination of the 
coefficient on the relative displacement of the water 
at the lake ends. This relative displacement is called 
the total wind tide, and is denoted by AH. The 
displacements are observed at Toledo and at Buffalo. 
In the case of a westerly gale if h B is the wind tide at 
Buffalo and h T is the wind tide at Toledo, the total 
wind tide is AH=h B — h T . 

This method of total wind tides has another advan- 
tage in the case of Lake Erie which may be mentioned 
very briefly. The particular phases of the storm 
effects that are specifically useful in determining the 
coefficient of wind stress are the maxima of the wind 
tides. Now, because the lake ends differ in shape, 
the rates of rise and fall of the water levels are differ- 
ent for the two ends. For Buffalo, the top part of the 
curve of tides is generally pointed; for Toledo, flat. 
Moreover, the maxima at Buffalo tend to occur 
earlier than at Toledo. These differences become 
more pronounced when the maturing storm increases 
rapidly in intensity. By adopting the total wind 
tide as the datum to be treated, the effects of these 
peculiarities of the individual behaviors are somewhat 
reduced and thus the data become more reliable. 

In some of the charts of the westerly gales received 
from the Lake Survey the displacements of water for 
the western end of the lake are shown for Gibraltar 
or Amherstburg. In these cases the wind tide indi- 
cations shown had to be reduced to the Toledo values. 
For the reductions it was necessary first to establish 
the correlation between the displacements at these 
western localities. This was possible, since in some 
of the charts the simultaneous values of the displace- 
ments are shown. Also Blunt [2] gives values of the 
maximum displacements for these localities. Intro- 
ducing all these in figure 9 and drawing a straight 
line through the plotted points, it is seen that on the 
average 

h T =1.71h A (12) 

where h A is the displacement of water either at 
Gibraltar or Amherstburg. This is the correlation 
to be applied. 

A final matter to be considered relative to the 
observed values of the wind tides is the correction for 
the effect of barometric pressures. Ordinarily, when 
the gale is moving over the lake area the pressures at 
the two ends of the lake are unequal. When the gale 
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Figure 9. 



O, Oct. 14, 1893. 
• Nov. 22, 1893. 
D, Nov. 28, 1893. 
B, Oct. 28, 1895. 
O, Feb. 11, 1896. 



♦ , Jan. 20, 1907. 
V, Sept. 25, 1941. 
T, Jan. 2, 1942. 
>, Jan. 3, 1942. 
► . Mar. 25, 1947. 



Data for years 1893, 1895, 1896 from Blunt. 

is a westerly, the greater pressure is at Buffalo. The 
difference of the atmospheric pressure at Toledo and 
Buffalo every hour of the storm period may be 
expressed as a head of water to be subtracted from 
the observed total wind tide. 

5. Data for Winds and Wind Tides of Lake 
Erie Storms 

The curves of wind intensity and total wind tide 
for the storm of December 31, 1911, shown in figure 8 
are typical of all the storms considered in the present- 
study. Taking the presentation in figure 8 as an| 
example, it may be seen that the entire manifestation 
of the storm may be broken into three epochs; the 
epoch of maturing storm, the epoch of relative 
steadiness, and the epoch of recession. The hydro- 
dynamical behavior of the w^ater in the lake is 
expected to show marked differences during thesr 
epochs. In the initial epoch the wind must blow for 
some length of time before the response of the watei 
to the action of the wind is completed. The reason: 
for this condition is that any manifestation of wind 
tide must be associated with the flow of water from 
one end of the lake to the other. The action of the 
wind must establish a layer of drift current, the 
thickness of which must increase with time either 
under action of molecular viscosity, or turbulent 
viscosity, or both. As the water is being collected 
at the leeward end of the lake, a returning gravity 
current is created which likewise increases in intensity 
with time. The second epoch represents that steady 
condition in which the transport of water through 
the body of the drift current is counterbalanced by 
the transport of water in the returning gravity current 
that is maintained by the unchanging surface 
gradient of the lake waters. The third epoch is 
associated with decreasing wind intensity. With the 
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strength of the wind decreasing, the stress of the 
\ id cannot maintain the adverse gradient of the 
iurfaee water with the result that a surge of water in 
bvave motion with a period equal to the seiche period 
of the lake is produced, in some records these 
seiches with periods of 4 hr are clearly seen during 
the time of the falling storm. The inertia effects are 
very pronounced in the extreme epochs: in the 
intermediate epoch they are considerably reduced. 

For the determination of the relation between 
wind velocity and wind stress, only the manifestation 
of water level changes in the intermediate epoch will 
be considered. Referring once again to figure 8, 
the reference time is the instant of maximum wind 
tide. In each case of the storms considered the 
reference time is evaluated by taking the average of 
I the times of maximum deflection of water at the two 
ends of the lake. As the total wind tide of the 
intermediate epoch we have taken the average value 
!of the wind tide over a duration of 4 hr around the 
reference time. For the corresponding wind values 
we have taken the average value over a dn ration of 
5 hr immediately preceding the reference time. 

The wind velocities and wind tides thus obtained 
from each record are entered in table 2 with the dates 
of the storms indicated. Tn some of the wind tide 
"harts the storms showed two peaks with an extended 
flat region between them. The total tides of such 
regions are also included in the data of table 2, and 
are identified with asterisks. 

Table 2. Wind tides of Lake Erie during westerly gales 



No. 


Date 


V 


V-' 


Ml 






in ph 


(mphY- 


ft 


1 


Nov. 21, L900 


50.5 


2549 


13. 12 


2 


Oct. 20. L905 


31. l 


983 


6. os 


3 


Oct. 2(). 1906 


38. 3 


1471 


9.75 


4 


Jan. 20. 1907 


is. 1 


2312 


12.04 


5 


Dec. 7. looo 


40.1 


ions 


10.51 


6 


Dec. 31,1911 


38. 9 


1516 


9. 53 


7 


Jan. 31,1914 


34.7 


i207 


7.95 


8 


Dec. 0.1017 


43.2 


1865 


10.17 


9* 


Dec. o.ioiT 


33. 1 


1093 


4.56 


10 


\)vc. 10,1907 


34.3 


1170 


7. 62 


11 


Dec. 18, hum 


45. 1 


2033 


L2.30 


12 


Dec. 8, 1027 


47.4 


2243 


L3.24 


13 


Dec. 9, 1027 


27.2 


712 


1. 13 


14* 


Dec. 9. 1927 


20. 1 


684 


3.44 


15* 


Dec. 9, 1927 


22.1 


487 


1.73 


16 


April 1,1929 


51.3 


2634 


L3.31 


17 


Jan. 241939 


38.8 


1508 


9.40 


18 


Sept. 25.1941 


35.7 


1277 


9.06 


19 


Jan. 2, 1942 


40.4 


1632 


12. 53 


20 


Jan. 3. 1042 


19.5 


379 


2.38 


:J 


Nov. 22, 1946 


34. 2 


1171 


8.36 


22 


Mar. 25, 1947 


35.7 


1276 


8.34 



The square of the wind velocities, V 2 , near the 
ake surface and the corresponding wind tides, AH, 
taken from table 2 are plotted in figure 10. Com- 
paring the position of the plotted points with respect 
to the full line drawn, it is seen that the data cor- 
responding to high wind velocities show small 
peviations from the curve, whereas those corre- 
sponding to low wind velocities show proportionately 
larger deviations, and all towards lower values of 
AH. One inference is that the wind stress for low 
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Figu RE 1 0. Wind tides in Lake Erie. 
Data for plotted points taken from table 2. 

wind velocity is relatively smaller. Unfortunately, 
the observations for the lower wind velocities are not 
sufficiently numerous to establish this point bevond 
doubt. 

6. Integration of the Wind Tide Equation for 
Lake Erie 

First the wind tide differential equation for a lake 
of arbitrary shape will be derived. Let positive x 
be drawn along the direction of the wind, this direc- 
tion coinciding with the axis of the lake. At the 
point x the wind tide is h, it being the elevation of 
the raised water surface with respect to the level of 
the undisturbed water. See figure 11. Corres- 
ponding to h at x the surface width of the water is 
B, the cross-sectional area is A, and the mean depth 
is H. The cross section is normal to the axis of the 
lake. The velocity of the wind is V, and the cor- 
responding wind stress is r s . At the bottom the 
stress induced by the gravity current is r . The 
density of water is p and that of air is p a . Two sec- 
tions Ax apart, one at x and the other at x-\-Ax, are 
taken. If the total pressure force on the section at 
x be denoted by P, the pressure force on the section 
at#+A# will be P+Apg (dh/dx)Ax. Hence, the dif- 
ference of the two pressure forces is Apg (dhjdx)Ax. 
The inertial action of flow in the drift current above 
and in the gravity current below being ignored, the 




Figure 11. Notation diagram for wind tides. 
II is the mean depth in the cross section of the lake at x 
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only remaining counterbalancing forces are the frac- 
tional forces amount to Z?(t s + t )A2?. Thus, equat- 
ing the forces and dividing by Ax, Apg (dh/dx)=B 
(t s +t q ), and since A = B (H+h), then 



dh 



Ts+Tq 



dx pg(H+h)' 



(13) 



which is the wind tide equation for a lake of arbi- 
trary shape. As the volume of water in the lake is 
not changed, the condition of continuity reduces to 



Bhdx = 0, 



(14) 



where L is the length of the lake axis. 

Equations 13 and 14 may be changed into dimen- 
sionless forms. The average section depth may be 
expressed in terms of the average lake depth H , 



ft-J-f 

BLJo 



L BHdx. 



(15) 



The section surface width may be expressed in terms 
of the average lake width, B, 



B=~ f L Bdx. 



(16) 



The distances along the axis may be expressed in 
terms of the lake length, L. 

Putting 



and 

eq. 13 and 14 reduce to 

d 1= _H /H 

di (l 
and 



T s + T L 

pgH H ' 


(17) 


h 
V= H Q ' 


(18) 


i-- x 
»-7.' 


(19) 



<n (l+uflo/ff) 



;§•*■ 



=0. 



(20) 



(21) 



It is important to remark that a is the parameter of 
stress, rj, the parameter of wind tide, f , the parameter 
of distance. In a lake of arbitrary shape H /H and 
BjB are functions of f . 

As ordinarily rj is a small quantity, eq 20 may be 
written in the form 

d v H / H \ 

ri^iv (22) 



The solution of the equation may be expressed in the 
form 

(23) 
(24) 



where 



Jo H 



(25) 



the quantities C x and C 2 being the constants of in- 
tegration. Now, from eq 21 



and 



r 1 B 

Jo B 

r^F 2 u)d^=o, 

Jo B 



(26) 



(27y 



and these serve to evaluate the constants C\ and C 2 . 
To carry out the determinations of the functions 
Fi and F 2 , it is necessary to obtain the shape quanti- 
ties H/H and B/B for Lake Erie. These quantities 
as determined from the Lake Survey Office map of* 
Lake Erie are shown in figures 12 and 13. Using the 
data from the figures, and following the method of 
computation shown above, the functions F x and F 2 
are determined. These are plotted in figure 14. For 
the subsequent analysis the end values 



^(0) = -0.641, #>(0) = 0.323, 

F 1 (l)= 0.512, ^(1) =0.118, 
are important. 



(28) 
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Figure 12. 
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Mean depths of Lake Erie as function of distance 
from Toledo, 
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Figure 13. Surface widths of Lake Erie as function of distance 
from Toledo. 
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Figure 14. Values of the functions F\ and F> for Lake Erie. 



At this point, some observed values of the lake- 
water displacement may be compared with the theo- 
retically determined curve of the analysis above. 
The most suitable way to represent the curve is to 
consider the ratio 2h/(hi—h ), where h is the surface 
displacement at x, and h and hi are the extreme 
values, as a function of x/L. Since c 2 is a small 
quantity, F 2 may be neglected and thus, 



2h 



2*\(f) 



K-K F.W-F^O) 



(29) 
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Figure 1">. Comparison of computed wind tides for Lakt Erie 

with observations. 

a, Oct. 29, 1892 (Blunt). O, Oct. 14. 1893 (Blunt). T. Sept, 26, L941. • 
Jan. 2, 1942. ■, Jan. 3, 1942. 



The quantity on the right hand, computed with 
reference bo the values in figure 1.4, is shown in figure 
15 by the full curve. The plotted points are from 
observations of different storms. The points repre- 
sented by open circles and open squares are taken 
from data given by Blunt [2]. Apparently there 
appears to be a fairly good agreement \ bei ween 
theoretical and observed values. 



7. Evaluation of the Coefficient of Wind 
Stress 

Representing the end values of t\ as tjo and rj h it is 
seen from eq 21 that ^ =<r^\(0)— a 2 F 2 (0), and r} l = 
aF.il) -a 2 F 2 (l). Hence ( Vl -r }o )=Ar ) = a[F 1 (l)- 
F 1 (0)] — o- 2 [F 2 (l)—F 2 (0)] } and introducing the values 
in eq 28, Ai?= 1.153 <r + 0.205 a 2 , or o-=0.867 An— 
0.178 a 2 . Treating the second term in the right hand 
member as a small quantity, 

cr=0.867A7 7 -0.134Ar. 

Introducing values of <j and n from eq 17 and eq 18, 
and remembering that AH=h\ — ho, we have 



T s +Tq L _ AH 

pgH H Q ' 0q 



—(f)' 



(30) 



This is the relation that relates the wind stress to the 
total wind tide for Lake Erie. 
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Expressing the shear at 
fraction of the wind stress, 



the lake bottom as a 



To = UT S) 



(31) 



the relation in eq 30 may be written as 



1+n 



[0.867 f -0.134 (fJ]^ p(J H .(32) 



Introducing the definition of the coefficient of wind 
stress, 

T s =XPaV 2 , (33) 

where p a is the density of air, and V the velocity of 
the wind, we have from eq 32 

x-i+nl 1 ° 16 H JTT 7al^' (34) 

which is the proper formula to evaluate the coefficient 
of wind stress from the total wind tide observed at 
Lake Erie. In the application of the formula the 
depth of water H is taken as 58 ft., the ratio of the 
density of water to the density of air, p/p a , as 1,250. 
The wind velocities V are to be expressed in feet per 
second and the wind tides AH in feet. It may be 
supposed that the frictional stress of the lake bottom 
is about one-tenth of the wind stress. Accordingly, 
n is 0.10. 

Applying the above formula to the set of wind 
intensities and corresponding wind tides appearing 
in table 1, a set of values of the wind stress coefficient 
are obtained which are entered in table 3. The same 

Table 3. Coefficient of wind stress deduced from Lake Erie 

gales 



No. 


V 


x x 103 


No. 


V 


x x 103 




mph 






mph 




1 


50.5 


2. 13 


12 


47.4 


2.44 


2 


31.4 


2.87 


13 


27.2 


2.35 


3 


38. 3 


2.78 


14 


26.1 


2.15 


4 


48.1 


2. 17 


15 


22. 1 


1. 53 


5 


40.1 


2.73 


16 


51.3 


2.88 


6 


38.9 


2.63 


17 


38.8 


2.64 


7 


34.7 


2.76 


18 


35.7 


2.96 


8 


43.2 


2.23 


19 


40.4 


3.22 


9 


33. 1 


1.78 


20 


19.5 


2.68 


10 


34.3 


2.71 


21 


34.2 


2.99 


11 


45.1 


2.52 


22 


35.7 


2.77 



data are plotted in figure 16 as full circles. Plotted 
in the same figure as open circles are the wind stress 
coefficients which Neumann obtained in an exami- 
nation of the Pelman observations for the Gulf of 
Bothnia [7]. Neumann's original values were reduced 
to correspond to winds 25 ft above the surface of the 
water. 

An examination of figure 16 reveals that in the 
region of the higher wind velocities the coefficient of 
wind stress is practically independent of the velocity 
of the wind. Furthermore, the coefficient of wind stress 
for the two bodies of water are of like value. The 
average value of the individual determinations is 
0.00254 for Lake Erie, and 0.00236 for the Gulf of 
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Figure 16. Comparison of coefficients of wind stress in Lake 
Erie and the Gulf of Bothnia. 
• . Lake Erie. O, Gulf of Bothnia. 

Bothnia. Now, the data from the Gulf of Bothnia 
for low wind velocities indicate that for this region 
the coefficient of wind stress is not constant but 
decreases with increasing wind. At the moment no 
similar data are at our disposal for Lake Erie and 
efforts are being made to remedy this gap in our 
investigation. 

8. Evaluation of Sea Roughness 

We shall adopt for the determination of sea 
roughness the procedure used by Neumann. Writing 
eq 33 in the form (V/V*) 2 =x~ 1 , substituting in 
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Figure 17. Variation of sea roughness, €1, with wind ve- 
locity, V 
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eq 1 f and changing e to ei, there is obtained 
-==5.75 lo£>* :r -^- 



(35) 



which relates the sea roughness e x to the elevation 
z n for which the wind is determined and the co- 
rd icient of wind stress x- In the present case, since 
X is determined from wind velocities prevailing at 
a height of 25 ft, 2 n = 25 ft, or 762 cm. As mentioned 
a hove, the average value of the coefficient of wind 
stress from the Lake Erie data is 0.00254. This 
yields for the sea roughness the value e^O.27 cm. 
As will be remembered, in the computations required 
to reduce the city wind velocities to the lake velocities 
it was necessary to assume the value of the sea 
roughness beforehand. The assumed value was 
€i = 0.3 cm, and this selection is now seen to be 
satisfactory. 

The full line drawn through the data points of 
figure 16 may now be used to investigate the variation 
of sea roughness with wind velocity. Using values 
from the curve and the formula above, eq 35, the 
results of the computations are given in figure 17 
in the form of a curve. With increasing wind 
velocities the magnitude of the roughness decreases. 
When the higher wind velocities are attained, the 



roughness approaches a constant value, pract ically. 

We refrain at the present moment from discussing 
the theoretical significance* of the facts that both the 
coefficient of wind stress and the sea roughness 
decrease with increasing wind velocities. 
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